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ELECTRON TUNNELING AS A BASIS FOR SEMICONDUCTION IN PROTEINS 
by Lawrence Flax and Dennis Flood 
Lewis Research Center 

SUMMARY 

Organic solids composed of biological macromolecules such as proteins are known 
to be semiconductors. Electron tunneling from molecule to molecule in the solid has 
been proposed by others as a possible mechanism for the electrical conduction which is 
observed. Apparent qualitative agreement between the form of the experimental and 
theoretical conductivities based on tunneling has been shown; but no attempt to obtain a 
quantitative estimate has, until now, been made. This report is a quantitative investi- 
gation of such a mechanism within the framework of the WKB approximation. The inter- 
molecular potential barrier is represented by the sum of the image charge potential and 
an applied potential which produces a bias between neighboring molecules. An approx- 
imate form for the tunneling coefficient is developed which is valid for applied electric 
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fields below about 10 volts per centimeter, a limit well above normal experimental 
conditions. The conductivity, calculated on the basis of the model and assumptions 
chosen, is several orders of magnitude smaller than that observed experimentally. The 
severe quantitative disagreement arises from the fact that the potential barrier sym- 
metry is relatively unchanged in low electric fields. In addition to the quantitative fail- 
ure, the tunneling model also fails to predict the well-known T Compensation law" for 
the behavior of the conductivity with activation energy. These results suggest that elec- 
tron tunneling per se is not the mechanism for the conduction of electrical current in 
proteins. 


INTRODUCTION 

Semiconduction in materials composed of organic macromolecules is an established 
empirical fact. There exists, however, no suitable theoretical explanation which can 
account for many of the experimental details that have been observed. For example, 
while the temperature dependence of the conductivity is given by 



-E A /2kT 


( 1 ) 


ct(T) = <j Q e 

where the activation energy E A is in the range of 2 to 3 electron volts, the so- 
called "pre- exponential factor, " a Q actually has the form 

a c = e^ E A + ^ (2) 

where £, and t] are independent of both T and E A - Equation (2) is known as the com- 
pensation law. Eley (ref. 1) suggested some possible mechanisms which might explain 
the behavior of the conductivity in this form, but did not investigate them extensively. 
Following one of the suggestions, Kemeny and Rosenberg (ref. 2) treated an organic 
molecular solid as a collection of molecules with potential barriers surrounding them, 
and calculated the conductivity by assuming that the electrons tunnel from one molecule 
to the next. They considered only the simplified cases of either rectangular or triangu- 
lar intermolecular potential barriers, and obtained qualitative agreement with the com- 
bined form of equations (1) and (2). 

This report presents a calculation of the tunneling current using an intermolecular 
potential barrier which we believe is more realistic. An expression for the current 
density is developed within the independent particle and WKB approximations, and a 
numerical estimate of the magnitude of the current density is made. 


ELECTRON TUNNELING 

The probability per unit time that an electron in state a on one side of a barrier 
will make a transition to state f3 on the other side is (ref. 3) 

<3) 

where is the density of states for one spin at /3, and f ff and fg are the probabili- 
ties of occupation of the states a and /3, respectively. The quantity is the 

tunneling matrix element for a transition from state a to /3. Bardeen (ref. 4) has 
shown that the matrix element may be written as 

M a0 = - < 4 > 

where J Q( ^(x) is the matrix element of the x- component of the current density operator 
J evaluated in the barrier region: 
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•V*) = <f'J J xlV 


(5) 


where \{/ a and ^ are the wave functions which describe the states a and /3. The 
matrix element vanishes imless the transverse wave number kj. is conserved in 

the process. If state a is considered to be on the left of the barrier, the current from 
left to right is obtained by summing over all states a of fixed k^, summing over 

k^, multiplying by the electronic charge e, and finally multiplying by 2 to account for 
electron spin. Since the probability per unit time for a transition from state j3 to state 
a is 



(6) 


the net current density flowing from left to right through the barrier is 

J = 2 e Z/ (I V’« - w dE / Kfii V « (1 ■ 

k t h 

■ iMualV^ 1 -'«)<>»}“ < 7 > 

where E is the total energy for a fixed kj.. The density- of- states factors allow the 
sum over the initial states to be replaced by an integral. Furthermore, 


so equation (7) reduces to 


J ^ J" p a p j3^a " f /3 ^ dE 
k t 

The matrix element has been evaluated by Harrison (ref. 3) in the WKB approximation 
and is given by 

M “0 " (£) P a /2p f? 1/2 exp 2 M *) <10) 

where x^ and X2 are the classical turning points of the motion (see fig. 1), 
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( 11 ) 


k x = j^T [ v(x) - E x]} 1/2 

V(x) is the barrier potential energy, and E is the "x-directed” part of the electronic 
energy of motion. 


DESCRIPTION OF MODEL 


A sketch of the potential barrier employed in this calculation is shown in figure 1. 


V 



Figure 1. - Assumed electronic potential energy as a func- 
tion of distance within the organic solid. Protein molecules 
are represented by a series of potential wells surrounded by 
finite potential barriers. The bottom of the well is the 
electronic ground state. A band of unfilled, higher lying 
energy states is shown for illustration. Scales of both 
ordinate and abscissa are arbitrary. 


The organic solid is assumed to consist of a collection of molecules, each represented 
by a potential well surrounded by a potential barrier. Within an individual molecule the 
potential energy of an electron is assumed, as a first approximation, to be 

V(x) - 0 (12) 

The nature of the potential energy within the molecules does not affect the tunneling 
coefficient in any way, so that this assumption is unimportant at this point. (The poten- 
tial energy is important, however, in determining the flux of electrons impinging on the 
barrier. ) The space between molecules is divided into two regions. In the first region 
0 < x ^ (a/2), and in the second (a/2) < x < a, where a is the intermolecular separa- 
tion. In the first region the potential barrier may be approximated by 

V(x) - E t - eFx - 2 (13) 

x 
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and in the second 


V(x) = E t - eFx y— (14) 

a - x 

F is the externally applied electric field, Erp is the total depth of the potential well, 
and y represents the strength of the interaction between an electron and the charge 
induced locally on the molecule when the electron emerges from it. In each region, 
what has been done is to ignore the attractive force exerted on an electron by any mole- 
cule which is more than one- half the interm ole cular distance away. The forces cancel 
exactly when the electron is midway between the molecules (assuming that the molecular 
orientation is such that y is the same for any two neighboring molecules). Otherwise 
the forces tend to lower the barrier a small amount at points away from the midpoint. 


CALCULATION OF MATRIX ELEMENT 

The matrix element describing the transition of an electron from one molecule to 
the next has the form, in this case, 



The classical turning points are found by setting the integrands equal to zero. Thus for 
the region 0 ^ x s (a/2) 


E t - E x - eFx - = 0 


(16) 


for which the smaller root is 


E T~ E x 


X 1 = 


- |/ (E t - E x ) 2 - 4yeF 


2eF 


( 17 ) 
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Similarly, for (a/2) < x ^ a 


(Erp - E x )(a - x) 

- eFx(a - x) - y = 0 

(18) 

from which 



E t - E x + eFa - 

xr — . .. 

y (Erp - E x + eFa) + 4yeF 

(19) 

x 2 ~ 

2eF 

It will be convenient to make the following substitutions and change of variable 


e t - 

■ E x =b 

(20a) 

eF = c 

(20b) 

x=t 2 , 

0 s x ^ (a/2) 

(20c) 

a - x = u 2 , 

(a/2) s x < a 

(20d) 

The first integral thus becomes 



j 

*1 

J- t 4 + £ t 2 - 1 dt 

f c C 

(21) 

while the second is 



r C U 2 
I 2 =-2^ / 

•W2 

</u^ + - ca ) u 2 - — du 

fee 

(22) 

where 



t l : 


(23a) 

u 2 = 1 

l /a - x 2 

(23b) 


The integral in equation (21) may be cast into standard elliptic form by making the 
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substitutions 



Comparing equation (24a) with equation (17) reveals that 

Xj = p 2 


so that equation (21) becomes 


4 = 2 Vc f* 7 * V (t 2 - 


p 2 )(q 2 - t 2 ) dt 


(24a) 


(24b) 


(25) 


(26) 


This integral can be found in a table of elliptic integrals (ref. 5) and may be written as 

r 




2\/c 


Mq 2 + p 2 )E(<?,k) - 2p 2 F(<p,k)] 


a „2 2 

— q _ p 

2 


Vf 


vf 'IM 


> (27) 


2 2 

provided that q ^ (a/2) s. p . E(<p,k) and F(<p,k) are the incomplete elliptic integrals 
of the first and second kinds, respectively, and 
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By making the substitutions 


r 


2 


s 


2 



in equation (22) and noting that 


u b =a- 


x 2 


= s 


2 


I 2 becomes 


I 2 = -2VH 



+ s 2 )(u 2 


+ r 2 ) du 


(29a) 


(29b) 


(30) 


(31) 


This too is in a standard form and can be found in the tables (ref. 5). The result is 



(32) 


2 2 

provided that r 21 (a/2) s. s , and where 




r + s 


(33a) 


cos (p i = — — (33b) 

It should be noted that the two elliptic integrals are not quite the same. This arises 
because of the change in the origin of the dominant image charge as the electron tra- 
verses the barrier region. 
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Before inserting these results into the expression for it will be helpful to 

define two dimensionless variables, y and y^. They are 


y = 



*1 = 



(b - ac)^_ 


(34a) 

(34b) 


The matrix element may now be written as 

v(y) + ( l - ■?) Vv w(yi) ]} <35) 

where 


M 


a/3 


= ± n~ 1//2 n-J/ 2 exp - JV ^ 


2 77 


a 


3fieF 


(Er 


V 


3/2 



and 
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CALCULATION OF CURRENT DENSITY 

The expression for the current density is, after inserting equation (35) into equa- 
tions (9) and (15) 



where Ej is the energy associated with the motion of an electron parallel to the barrier 
surface (i.e. , perpendicular to the net current direction) and is a function of k^.. The 
basic assumption proposed by Kemeny and Rosenberg (ref. 2) was that electrons are 
thermally excited from their molecular ground state to some intermediate state below 
E,p from which they then tunnel. Accordingly, the occupation probability of such a 
state on the left of the barrier can be assumed to be 



and on the right side it is 


Fp = exp 


/_ E + eFa\ 
\ kT / 


(40) 


Hence 


J 




E + eFa\ 
kT / 


dE 


(41) 


where D(E^, - E x ) contains all the other factors in equation (38) and is known as the 
tunneling coefficient. Converting the sum over to an integral and separating the 
"perpendicular" energy from the "x-directed" energy fields 
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J = 


m o 


°( E T " E X> exp ( 


2 2 ' 
‘ p y - p z 


2mkT 


exp 


•E 


kT / 


- exp - 


E x + eFa ^ 
kT j 


dE x d Py dp z (42) 


Carrying out the integration over p^ and p z , the momentum components in the y and 
z directions, respectively, 


47rmekT 


/ 

t/0 


D(Er 


E x ) exp | 


H 

- exp (- eFa )l 

\ kT/ 

V kT / 


dE. 


(43) 


The applied electric fields are typically less than 1000 volts per centimeter (~3 stat V/ 
cm). If y is assumed to be e^/4 (its value for the image charge induced by an electron 
leaving a metallic surface) and if Erj, - E x is of the order of 1 or 2 electron volts, 


y 2 =■ 


— ? as 5X10“ 6 F 


(44) 


(B t - E x )^ 


(F in stat V/cm). A representative range for the inter molecular separation a is 10 
7 

to 10 centimeter so that 


,-8 


y l = " 


e 3 F 


a* 5X10' F 


(45) 


(E t - E x - eFa)" 


as well. For values of y and y^ this small. 

v(y) ^ w(y 1 ) ^ l 

so that 


(46) 


D(E t - E x ) as exp 


- 4y2m / F F \3/2 

Z \F-j rp “ Hi ) 

3ReF 1 x 


(47) 
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The current density thus becomes 



The integration breaks into two ranges because D(E T - E x ) becomes unity when 
E 5: Erj,. The second integral is much smaller than the first and can be ignored. In 
addition, since eFa/kT « 1 for the range of field values chosen, equation (48) becomes 



Activation energies for semiconduction in many proteins have been found to be remark- 
ably uniform, falling within a range of 2. 6 to 3. 1 electron volts (ref. 6). This is con- 
siderably less than the ionization potentials of such molecules (E T ^ 10 eV). Instead of 
restricting the electrons to a single excited level, however, as did Kemeny and Rosen- 
berg (ref. 2), we assume here that a small band of energies exists, the lowest level of 
which lies in the observed range of activation energies. In such a case the logarithm of 
the tunneling coefficient may be expanded about the lowest level in a Taylor series. De- 
noting the lowest level by E Q 


D(E T ■ E x> K ■ E o > 3/2 - f < E T - E o> 1/2 < E x - V 


+ -< E l- E o>‘ 1 / 2 < E x- E o > 2 + 


3 

8 


(50) 


Inserting this expression in equation (49), 


12 



The integration in equation (51) should, in principle, extend only over the allowed band 
of energy- levels beginning at E Q . Because the integrand falls off so rapidly for higher 
energies, the error introduced by extending the range to infinity will not be too great. 

It will, if anything, produce an overestimate of the magnitude of the tunneling current. 
The integral in equation (51) may be readily evaluated (ref. 7). The result is 



13 


Considerable simplification of the preceding expression developed for the current 
density can be achieved by noting the relative magnitudes of some of the terms. A 
reasonable assumption for cp Q is 5 electron volts, while barrier heights are of the 
order of 6 to 10 electron volts, and E Q is of the order of 2 to 3 electron volts. Using 
the free electron value for the mass of the electron yields 


j/ 2m ?o _ 2.3X10 7 
KeF F 

At room temperature, 

P— « 1.2X10 13 
2kT 


Furthermore, 


8m< ?o _ 7.5X10 34 
3 (KeF) 2 F 2 


(54) 


(55) 


(56) 


so equation (53) becomes 


J 


47rme 2 Fa 

•y/2m 


, 1/2 


exp 


P 


exp 




3KeF 


KeF 


(57) 


j 


Since the argument of the error function is very large (~ 10°) the following asymptotic 
expansion may be used (ref. 7). 


(erf)(x) = 1 - (a + a g£ 2 + a^P) e ” X + 


(58) 


where 


1 

1 + nx 

(59a) 

n = 0. 4705 

(59b) 

a x = 0.3480 

(59c) 
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i 2 = -0. 0959 

(59d) 

a 3 = 0. 7479 

(59e) 


For an electric field of 300 volts per centimeter (1 stat V/cm) and values of electron 
mass and barrier height used previously, l ^ 3x10" . This means that 


J^T . x 

\ ReF J 




(60) 


Hence, the current density becomes 


J = 3 



e^F 2 a exp 



exp 




3ReF 


(61) 


from which the electrical conductivity is 


o 


F 


3 



e°Fa exp 




(62) 


DISCUSSION 

The conductivity given by equation (62) is not in the form required by the compensa- 
tion law. If, though, the assumption is made that 

3E \ 

— - + • • | (63 

2F t J 

then equation (62) can be rewritten to have the proper form. This is the approach used 
by Kemeny and Rosenberg (ref. 2) in their calculation of the conductivity using both 
rectangular and triangular interm ole cular potential barriers. Since E Q appears to be 
from 20 to 30 percent of E,p, however, stopping at the linear term is not obviously 
justified. Thus, there is a failure to agree even qualitatively with observed experi- 


^ /2 - <s T 


E o> 


3/2 


E^/ 2 
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mental behavior. In addition, the model fails to produce any degree of quantitative 
agreement. Again, for the relatively high applied field of 300 volts per centimeter, 

-log 10 a ** 1.3X10 6 (64) 

Observed values of -log^Q a are in the range from 0 to 20. The large value for the 
negative of the logarithm means that a is essentially zero. 

The failure of this model to predict the magnitude of the conductivity is not surpris- 
ing. For low values of the applied electric field (i.e. , <10 V/cm), the potential bar- 
rier will not be distorted much from its shape in zero applied field. Hence, the proba- 
bilities for tunneling from either side of the barrier to the other will be nearly equal, 
with the result that the net current will be extremely small. This fact does not appear 
to have been considered by Kemeny and Rosenberg in their discussion. The expressions 
they used are valid if tunneling is considered to proceed in one direction only. Such is 
the case, for example, in field emission from a metallic surface into vacuum, where 
there are filled electron states on only one side of the barrier. Tunneling through poten- 
tial barriers between neighboring molecules does not have this natural asymmetry in the 
direction of current flow, since presumably there are filled electron states on both sides 
of the barrier. It is possible that the states are not identical on each side of the barrier, 
which could introduce the necessary asymmetry, but that is a question which is beyond 
the scope of the present discussion. 

The restriction to low fields was made in this calculation only because the"? are nor- 
mally encountered in practice, and because considerable simplification of the calculation 
is achieved by so doing. An investigation of the behavior of the tunneling current at 
higher electric fields presently can only be done with the aid of a computer and is under 
investigation. The implicit assumption has been made that the electric field existing 
between neighboring molecules is just the applied field. The presence of some kind of 
field- enhancing effect which might alter the local electric fields could drastically affect 
the barrier shape, and therefore the current density. This possibility will also be ex- 
plored. 


CONCLUSIONS 

The conductivity obtained from the model used can be written in the form of the 
"compensation law, " but the approximation required to do so (eq. (63)) is not a particu- 
larly good one. More serious, however, is the complete lack of quantitative agreement 
with observed conductivities. This failure is directly attributable to the fact that the 
barrier symmetry is not altered appreciably by low electric fields, such as are normally 
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encountered experimentally, so that the net current through the barrier is essentially 
zero. The present results suggest, therefore, that electron tunneling, unassisted by 
any enhancing effects, is not a suitable mechanism for the explanation of the electrical 
conductivities of proteins. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 30, 1971, 

112 - 02 . 
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APPENDIX - SYMBOLS 


a 

intermolecular separation 

j 

X 

b 

difference between barrier 



height and x- directed 
energy, E T - E x 

J a/3< x > 

c 

product of electronic charge 



and applied electric field. 

k 


eF 

k, kj 

D(E t - E x ) 

tunneling coefficient 


E 

total energy of an electron 

k t 

e a 

activation energy 


E o 

electronic energy level to 
which tunneling occurs 

k x 

e t 

depth of molecular potential 

M a(3 


well 

m 

E x 

x-directed electronic energy 
of motion 

P o0 

E 1 

electronic energy of motion 


parallel to barrier surface 

p y 

E(<p,k) 

incomplete elliptic integral 

p z 


of first kind 

T 

e 

electronic charge 

t 

erf 

error function 

u 

F 

applied electric field 

V(x) 

F(<P,k) 

incomplete elliptic integral 



of second kind 

v(y), w(y 1 ) 

f *’ f /3 

occupation probabilities 

x l» x 2 

h 

Planck's constant 


R 

Planck's constant divided 

Y,Yl 


by 2n 

a, /3 

i 

V-i 

y 

J 

current density 



x- component of current 
density operator 

matrix element of 

x- component of current 
density operator 

Boltzmann constant 

moduli of incomplete elliptic 
integrals 

transverse component of 
wave number 

x- component of wave num- 
ber 

tunneling matrix element 

mass of an electron 

probability per unit time for 
transition from one elec- 
tronic state to another 

y- component of momentum 

z- component of momentum 

temperature 

integration variable, 

integration variable, -y/a - x 

intermolecular barrier po- 
tential energy 

dimensionless variables 

classical turning points of 
electronic motion 

dimensionless variables 

electronic state indices 


\ 

k 


18 


constant in image charge po- 
tential energy expression 



| , ?] constants 

P a , Pfi density- of- states factors 

ct(T) electrical conductivity as a 

function of temperature 

a Q pre- exponential factor in elec- 

trical conductivity 


<p Q difference between barrier 

height and energy level to 
which electrons tunnel, 

E T" E o 

<p , cp ^ arguments of incomplete ellip- 
tic integrals 




wave functions for electrons in 
states a and /3 
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